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.ABSTRACT 


In this thesis an attempt is made to model the signal 
received from a target, such as a submarine, by an active 
SODaR, such as a homing torpedo. For an iso-speed shallow 
water channel, the channel exhibits discrete multipath 
scattering. For typical bODAR beam patterns and normally 
encountered boundary reflection losses, the number of paths 
with significant energies contributing to the received 
signal is small. 

A characteristic feature of the above channel is seen 
to be the common doppler shift in the composite received 
signal due to target motion. This doppler shift is taken 
to be the average of all the doppler shifts of contributing 
rays as doppler ratio differences among paths are seen to 
be small. 

Modelling the attenuation coefficients as Zeromean 
complex Gaussian random variables and modelling the random 
travel times with appropriate density functions, the 
optimum receiver is found to be too complex to be of practical 
interest. If perfect knowledge of arrival times is assumed, 
the optimum receiver is of simple nature. A conventional 
receiver under this characterization suffers no marked degra- 
dation. With an appropriate choice of envelope of the trans- 
mitting pulse, the optimum receiver implementation is simpler. 
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CHAPTER 1 
INTRODUCTION 

The normally used receiver for the detection of a 
target in presence of white noise by an active SODAR (S und 
Detection And Ranging) is a correlation or a matched filter 
receiver. For such receivers the channel is considered to 
be dispersionless and unbounded, and the target is modelled 
as a slowly fluctuating point target. Then the return signal 
in an active SODAR is a replica of the transmitted pulse 
appropriately weighted and time delayed. Even if coloured 
noise is present, the same techniques used above suffice and 
it is generally observed that a suitable choice of trans- 
mitting waveforms will improve the receiver performance. 

? For bounded channels, there was no special mention 

on the performance changes for the above receivers. In this 
thesis we attempt to restructure the target and channel 
models when the channel is bounded and suggest some suitable 
receiver structures. 

A bounded channel is assumed to have well defined 
closely spaced boundaries, (an example could be the shallow 
water channel where the source is operating at high frequencies) 
and such boundaries are said to produce specular reflection. 
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at high frequencies it is appropriate to characterize the 
reflection process to be regular and deterministic. Now 
in such a channel when the source/receiver and target are at 
fixed locations, the acoustic transmission may be described 
by ray theory and the channel is said to have a discrete 
multipath structure. Hence the received signal could be 
a super-position of weighted and time delayed replices of 
the transmitted pulse arriving at receiver after the 
'reflection 1 at the target. In addition, a moving target 
will also introduce doppler shift on carrier frequency along 
each path. Hence the net effect of the channel and a target, 
which is modelled to have zeromean complex Gaussian reflection 
coefficient, is the spread of received energy over ambiguity 
plane. The amount of spreading of the received signal over 
ambiguity plane may be found by using ray tracing programmes 
even for a channel in which the sound transmission is a 
function of depth [l]. 

It is observed that the general interest in a multipath 
environment in a bounded channel is to investigate the 
coherence degradation between signals received at two 
remotely located sensors [2,3]. This scattering phenomena 
is some times applied in passive tracking of a moving 
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source in shallow waters [4,5]. An elegant model for the 
channel is suggested by Jacyna, Jacobson and Clark [5], 
where a passive source is tracked by two sensors. 

In this thesis we adopt the Jacyna, Jacobson and 
Clark's model for an active SODAR to detect a moving 
target. 

In certain occasions it is observed that special 
signal waveforms like coded pulse sequences or analog 
modulated waveforms (LFM etc.) are used in channels where 
thermal noise is a predominent component, for detection 
of targets. In this thesis we attempt to analyse the 
requirements of such waveforms when the channel is bounded. 

In Chapter 1 of this thesis the problem under study 
is proposed. In Chapter 2 , we first review various concepts 
that are involved for computing the parameters needed in a 
multipath received signal formulation, and adopt a model for 
detecting a moving target in a shallow water channel for an 
active oODAR. In Chapter 3, we attempt to study the optimum 
receiver for the received signal characterized by random 
attenuation coefficients and random time of arrivals. V»e 
then observe the simplicity in optimum receiver design 
having a knowledge of time of arrivals. For this situation 
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we consider a conventional receiver and study the amount 
of its degradation compared to the optimum receiver. The 
effect of choice of transmitting waveforms are studied. 

In the end we discuss the results. In Chapter 4, vie 
conclude the discussion with some suggestions for future 
work. 



CHAPTER 2 


MODELLING OF SHALLOW WATER TARGETS 

2,1 INTRODUCTION: 

In this chapter we are interested in developing 
a model for an under water SODAR target. Our concern is 
to characterize the signal received by an under water weapon 
system when it transmits a signal that is ’reflected’ from 
a moving target. Such a situation for example, is encoun~ 
tered when an active torpedo homes onto a submarine moving 
in continental shelfs characterized by shallow waters. 

The received signal features are dependent on -the 
channel parameters which modify the acoustic propagation of 
the transmitted signal. Boundary conditions and channel 
parameters may have complex influence on the acoustic 
propagation. In Sec. 2.2, we consider the effect of' various 
parameters $f the medium on wave propagation and the 
relevance of ray acoustics to the present problem. In Sec* 
2.3, we give expressions for travel time, spreading loss and 
reflection associated with an acoustic wave. In Sec. 2.4, 
we consider the multipath structure of the channel. The 
influence of a moving target and a general treatment to 
compute the return signal in a moving target environment is 
studied in Sec*. 2.5. Apd finally in Sec. 2*6 we dexive-ap 
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expression for the return signal for the above specified 
example and discuss the results obtained. Our account of 
the target model, in particular the influence of the source 
motion on the received acoustic field, closely follows the 
work of G.ivi. Jacyna, M.J. Jacobson and J.G. Clark [5]. 

2.2 UNDER-WATER ACOUSTIC PROPAGATION PRINCIPLES: 

The acoustic wave propagation in an underwater 
channel is influenced by regular and random inhomogeneities 
of the medium [ 6]v Refraction, duct propagation etc. are 
the phenomena due to regular inhomogeneities. The regular 
inhomogeneities refer to the spatial and temporal variations 
in the mean characteristics like temperature, salinity and 
pressure etc. in the medium. In ocean media these effects 
are considerable in the direction of depth and the ocean 
medium is often said to have a layered structure. The 
random inhomogeneities on the other hand, cause phenomena 
like scattering or reverbaration. Which of the above two 
inhomogeneities is dominant depends on the specific situation 
under study. In this thesis we .are mainly interested in 
regular inhomogeneities only. 

The regular variations in temperature, salinity and 
pressure influence the magnitude of the velocity of propa- 
gation of the acoustic wave and consequently the path 
traversed by the acoustic wave. 
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ye now briefly consider the basic principles which 
govern the acoustic propagation in the medium* 

The propagation of sound in an under water channel 
can be described mathematically by the solution of the wave 
equation using the appropriate boundary and initial condi- 
tions for a particular problem* In a homogeneous and 
quiescent medium the velocity potential 0 at a point on the 
wave obeys the wave equation given by [ 7 ] 

V 2 0 - K = 0 (2.1) 

c 6t 

where 1 c* is the velocity of propagation of the wave and 
c is given as 

c = (K/^ q ) 1/2 

where K is bulk modulus and is density of the medium. 

The velocity of the particle on the wave is given by 

v = V0 

and the pressure developed due to this velocity is given 
by 



where p Q is the density of the medium. 
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Under the assumption of small perturbations in an ambient 
state, the above wave equation can easily be seen to follow 
from: 

i ) th e principle of c o nserv at ion of mass: This principle 
states that in a fluid the net rate of increase of mass per 
unit volume is equal to the rate of change of density. If 
P(x,t) is the density of fluid at a point X and time t, 
v(x,t) is the fluid velocity at thatmpoint , the equation for 
conservation of mass is easily shown to be given by 

|f + V* (Pv) - 0 

ii ) The Euler’s equation of motion for a fluid : This equation 
states that the mass multiplied by accelaration of a fluid 
particle is equal to the net apparent forces exerted on it 

by it's environment and by external bodies. Neglecting the 

i 

body forces such as gravity per unit mass, the equation of 
motion is written as 

p • n = -Vp 

where p(x,t) is force per unit area or pressure. 

iii) The equation of sta te: It describes the relation 

between the pressure of a fluid, and its density and temperature. 
If an excess pressure developed due to an acoustic wave is 
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p, then this excess pressure is proportional to the fra- 
ctional change in density (Ap)* Hence we write 
Ap = K-A/° where K is the bulk modulus of the medium. 

There are two theoretical approaches to the solution 
of wave equation (2.1). One is called the normal mode 
theory [8 ], in which the wave propagation is described 
in terms of characteristic functions called normal modes, 
each of which is a solution of equation (2.1). In a bounded 
medium the standing wave pattern occurs due to super position 
of waves travelling in opposite directions towards boundaries. 
If H is the depth of the channel, then for all wavelength 
equal to where j = 1,2,3,..., a solution to wave 

equation is said to exist. In ocean media, the standing wave 
pattern occurs along the z— axis (depth) and the wave propa- 
gates along the x-axis (horizontal axis). The normal modes 
are combined additively to satisfy the boundary and initial 
conditions. With normal mode theory we obtain at a point 
and at a given time the combined effect of all wave paths. 
This analysis is relatively easy to handle when the number 
of modes is small. If the point at which the effect of the 
acoustic source is computed, is at a large distance (of the 
order of 50 Km) or the source frequency is small (less than 
1000 Hz) the resulting number of modes will be small and 
normal mode theory may be applied. 
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The other approach to solving equation (2.1) is by 
ray theory [9 ]. Ray theory represents the acoustic field 
at a point as a sum of ray contributions, each ray emanating 
from the source and reaching the point after a certain 
number of boundary reflections. In media in which the sound 
velocity changes appreciably over a small distance (say in 
the order a wavelength of the sound wave), ray theory is not 
applicable. Hence ray theory may be seen to be applicable 
to high frequencies or short wavelengths in a practical 
situation, hay tracing can be done using suitable computer 
programes for the energy loss, time of travel along each 
ray path [ 10 ], 

An underwater weapon generally has the following 
operating specif icationss 

i) Source frequency range : 10 - 100 KHz 

ii) Acquisition range : 2-3 KM 

iii) Operating channel depths : 200 - 400 meters. 

These factors allow us to use ray theory to analyse 
the received signal characteristics of an underwater 
weapon, 

2 . 3 RAY ACOUSTICS: 

We consider the wave equation in two dimensions, say 
x and z where x signifies the horizontal axis and z signifies 
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the depth, in a ocean medium. The positive direction of z 
is towards increasing depth. Let the acoustic source be 
producing a harmonic disturbance. :.'e consider a homogeneous 
medium in which velocity of propagation of the acoustic 
wave 1 c* is constant. Then the wave equation (2.1) is 
written as 



consider 

0(x,z,t) = y(z) exp j(+ax~tot) (2.3) 

to be a solution of equation (2.2). 

liquation (2.3) describes plane waves propagating in the 
direction of increasing x (+ sign) or decreasing x(- sign). 
The phase velocity of the wave is given by v and is 
equal to ~ . 


Substituting equation (2.3) in (2.2) gives the differential 
equation 



2 

ry 


o 


where 



a 


2 ,v 
a (~j 

c 


1 ) 


(2*4) 


(2.5) 
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Consider a solution of (2.4) as 

ny (z) = Ae^ z + Be”^ YZ (2.6) 

where A and 3 are constants of integration. 

Combining (2.6) and (2.3), we have 

0(x,z,t) = A e J<± ax+ r z -" t) + B e J(±«*~Yz-“t) (2.7) 

where the first term describes a wave travelling in the 

direction of increasing z, i.e., propagating down-wards 

and second term corresponds to a wave propagating upwards. 

2 

we write equation (2.5) as ^ =* a 2 + y 2 and define 

c 

k 2 4 (2.8) 

c 

where k is the wave number in the direction of propagation. 

It is the magnitude of the vector K whose direction is normal 
to the wave front i.e., K defines the ray direction. 

If we consider a plane wave described by an 
equation 

0(x,z,t) = A pCax+yz-ut) 

then this plane wave is moving in a direction K whose magnitude 
k is given by or inturn we can say that the plane wave 

is having a and y as its x and z components. 

Consider Fig. 2.1 and if 0 is the angle the direction 
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of plane wave makes with z-axis, then the horizontal component 
a = k sin© and vertical component y - k cosQ. Fox a homo- 
geneous medium we can see K is constant* That is the wave 
will , not alter it’s direction and it' s magnitude remains 
constant as it propagates. 

Now we derive a solution to the wave equation (2.2) in 
a medium where c is not a constant* 


For a source producing harmonic disturbance consider 
a solution to the wave equation (2.2) as 

0(x,z,t) » A(x,z) eJC^x.n-wt) (2.9) 

where S, a function of x and z to be determined, we substitute 
(2.9) into wave equation (2*2) and separate the real and the 
imaginary parts. Ignoring the time dependence, we have 


V^A - [(f§)^ + (§f)q A + K 2 A = 0 
a nd 


2 V A . V S + AV S = Q 

Rewriting (2.10) and (2.11) as 

2 , 


( H > 2 + c || )2 - * 2 


= 0 


and 


V 2 $ + § VA.75 = 0 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


(2.13) 
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then equation (2.12) implies 


f oS \ 2 


/ dS\ 2 

(*§y) 


K 


(2.14) 


Equation (2.14) is called Eikonal equation (corresponding to 
the wave equation) and S(x,z) must satisfy (2.14). S(x,z), 

O 

leading to a solution to the wave equation, if VA.V S is 

A 

small in equation (2.13) must also satisfy 


= 0 (2.15) 

The general form of the solution of equations (2.14) and 
(2.15) will be 


S(x,y) = c j, x + lf(K 2 ~ c^ 2 y y (2.16) 

where c^ is a constant to be determined. 

Equation (2.16) describes a plane wave front in which K is 
the wave number. It is the magnitude of a vector K whose 
direction describes the ray path. Let K make an angle 0 
with horizontal axis (x-axis), then we can show that 
c^ = K coscc. This implies that c^ is the direction cosine 
of K along x-axis. 

Having obtained S(x,z), we substitute equation (2,9) 
now into wave equation and obtain a second order differential 
equation in A(x,z). Solving for A(x,z), we obtain 0(x,z) 
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Conditions 


^ 2 


p . — — 0 and ^ ? a S is small are satisfied, if 
the fractional changes in medium parameters like velocity are 
small over a wavelength. These conditions are generally met 
when the source operates at high frequencies. 


The solution, (2.16) of eikonal equation is seen to be 
a plane wave with K 6 ~ • If velocity of sound ’ c* is a function 
of x and z, the gradient of the wave front and therefore of the 
ray path change with x and z. Now we derive a condition which 
the ray path has to satisfy at each point. 


Each point on the ray satisfies the eikonal equation 

(•If ) 2 + (||) 2 = K 2 = 4 ■ 

c 

An infinitesimal ray length can be described as 


do kdl = k(dx +> dz) 


1/2 


(2.17) 


using Fermet*s principle (which states that rays in physical 

space are paths of stationary time i.e., for small variations 

in path length between two fixed points, the resulting varia- 

P 2 

tions on the travel time integral / dt is zero); if the 
variation introduced on the path length between two fixed points 
we have S / do = 0, then we have from equation (2.17), for cd*0 

S / K dt =<?/ -f!= 0 (2.18) 

If we consider the velocity of wave 1 c' to be a function 

of z only, which is true generally in ocean media, then 
dx 

writing x z = ££ equation (2.18) implies 

2 

ctzT' 


y\ x 

s/ - z 


+ 1^ 


dz = 0 


(2.19) 
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Applying Euler-Lagranges formula, to equation (2.19) 


we have 


d .. _ dF 

r 

% z 


dir nr = 0 


( 2 . 20 ) 


where F 


V <* 2 + 1 } 

~cTzT" 


5F 


Equation (2.20) implies * = constant and we have 


1 

c 


z _ constant (=a) (2.21) 


f ( x z +1) 


From Figure 2.2, if 0 is the angle the ray makes with x-axis 
then 


= ~~ * cote 

Z GZ 


( 2 . 22 ) 


Therefore (2.21) implies, cos©/c( z) = a 

cos© _ 
cTz) - a 


(2.23) 


Equation (2.23) is called Snell's law and it is to be 
satisfied at each point on the ray. 


The eikonal equation in 3-u can be written as 

)y J 


(§) 2 + (f) 2 + (||) 2 - K 2 , 


If we assume the velocity of the wave ’ c ' is a function of 
one co-ordinate only, namely z -coordinate, then for s ray 
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path, the direction cosines OS/dx and dS/dy along x and y 
co-ordinates must be constant at all points. This implies 
f’b'S^^y} := consent and hence ray path is confined to a plane 
normal to xy plane, oo with no loss of generality we 
assume the ray plane to be xz plane. 


2,2.2 Computation of Travel Time and bpreading Loss: 


be have seen in a medium in which 1 c' is proportional 
to only z-co-ordinate, the rays are to be confined to xz 
plane only, be have also seen that at each point on the ray 
the Snell’ s law is satisfied, be now compute the horizontal 
disxance (along x-axis) a ray traverses when it travels from 
a given depth (say z = 0) to another depth (say z = z), 
Equation (2*22) implies dx = cot© dz, and the horizontal 
distance traversed is given as 

z 

x = / cot © dz (2.24) 

o 


p Z 

o f(l~z 2 c 2 (z)) 


dz 


(from equation (2.23)) 


If • c‘ is a linear function of z$ i.e, writing c = pz, where 
p is the gradient, then equation (2.24) is evaluated and 
we have 

X = const. +■ ( l-a 2 p 2 z 2 ) 1//2 (2.25) 
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Define II A and Xj ^ const,, the (2,25) implies 

(X - X x ) 2 + Z 2 = R 2 (2.26) 

Equation ( 2 . 26 ) implies that the acoustic ray traverses 

a circular path if T c' is a linear function of depth. The 

radius of the circular path is and the origin will lie 

ap 

at (X^,0) where z = 0 implies the depth where * c* is 
extrapolated to zero. Fig. 2.3 illustrates the above 
statement. 


Let 0 be the angle of inclination the ray makes at 

the source with the horizontal axis. Let c be the velocity 

o 

of propagation at the source. From equation (2.23) we have 


cos 0 

“cTzT 


a 


1 c 

Radius of curvature R = 

ap p cos© 

and the origin lies at a depth where velocity profile 

c 

extrapolates to zero or it is at a depth from the location 

X - '* 

of the source. Therefore centre of the circle is at (X.,0) 

c x 

where X. = tan©. 

1 P 

lie define travel time along the path is equal to 


~~ where dS = “v-rs as may seen from Fig. 2.2. 

Therefore if c is a linear function of z (i.e. , c = pz) then 
travel time 1 t* is given as 
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t = 


P 2 f(l-a 2 p 2 z 5 J) 


( 2 . 27 ) 


Transmission loss alone; a ray path is a sum of 
spreading loss and attenuation loss. Attenuation loss 
includes losses- due to absorption and scattering. In this 
discussion we ignore attenuation losses. Spreading loss is 
a geometrical effect representing the regular weakening of 
signal as it spreads outward from the source. In shallow 
water channels the medium is bounded by two surfaces and 
the spreading loss is assumed to be cylindrical loss and 
is given as [> ] 



where 

p is the source energy per unit time per unit solid 
angle . 

dA is the area swept oat by the wave surface, normal 
to the ray at a horizontal distance x q then 



x sin© 6x i-l 

n smrnarnm ifrWtaBxai I 

p cose o oe o 1 


( 2 . 28 ) 


whore 0 Q is the angle the ray path makes with the horizontal 
at the source and 9 corresponds to the angle the ray makes with 
the horizontal at the point of interest. 
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2.3,2 Reflection Coefficient: 

At the interface separating two media, the ray has 
to satisfy the following ooundary conditions: 

i) the vertical particle velocity must be continuous. 

c>0 d0 o 

= < 2 - 29 > 


ii) continuity of pressure across the boundary is to be 
maintained, i.e. 


p b ±\- p 
1 1 


^2 

Tt = * 2 ™3T 


•(2.30) 


where P ^ and are "the densities of the media. 


Define 0^ n = A exp (ax + yz - wt) as incident ray and 

0 = B exp (ax - yz - wt) as the reflected ray. Refer 

ye 

Fig. 2.4. 


Substituting above equations into boundary conditions, 
the reflection coefficient turns out to be 


B 

A 


/°2 

<*1 


~ ¥ 


( C l/ C 2 )2 “ Cos2e i7 

. ' ““ "■'”2 ! 

1 - cos 9, - 1 


Vv- 


<°1 


+ f 


(Ci/Co) 2 - cos 2 ©, i Va- 


'1 12 _ _ 
" r o 

1 - cos ©. 




(2.31) 
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where c^ and are the velocity of propagation and 
density of ocean medium and c^ and corresponds to same 
quantities in the second medium. 


Consider the second medium to be air and the interface 
to be ocean surface. Now the interface separates two media 
having large difference in their acoustic impedance. Acoustic 
impedance is defined as P c, and we have /° c^ = l.oxlcP g/ 
(cm) sec. and ~ gm/(cm; sec. (for air), 

bo equation (2.3l) implies that the reflection coefficient is 
unity and the reflected ray will have it’s phase advanced 
by 180°. 

On the other hand at the interface separating the 
sea and the sea bottom with c 2 as P ara meters corres- 

ponding to the sea bottom, we have c^ < c^ and P ^ < P^. 

yve define critical angle © as 

J c 

C, 

cos © c = (2.32) 


Now from equation (2.3l) and (2.32) we see that the 
reflection coefficient will increase from a value equal 
to for 9 = 90° to a value, unity at angle of 

P 2 c 2 + Pl c l 

inclination © = 0 . A typical curve for reflectivity vs. 

V 

inclination angle is shown in Fig. 2.5. For the angle of 
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inclination © from © to 0°, the reflection coefficient is 
purely imaginary with magnitude of the reflection coefficient 
being unity and the phase change is determined by inclination 
angle. The phase angle is given as 


6=2 tan 


-1 


f> 1 f cos 2 ©-(c 1 / c 0 ) 2 

K 


'/o 


r cos 0-(C./C 9 ) ~J 

/ — — -v*— (2 - 33 > 

L i - cos^e J 


and £ is 0° at the critical angle and increases to 180° 
at 0= 0°. 

2.3,3 Ray Displacement: 

In general, the fields computed by the ray theory 
differ from those computed by the normal mode theory. Tinble 
and Bold [l2j account for this discrepancy by the beam 
displacement suffered by rays on total reflection from the 
bottom. Total reflection means that for rays whose inclin- 
ation angle are less than critical angle, the reflection 
coefficient amplitude is unity and has only phase change. 

And this additional phase is considered as extra path. When 
the extra path lengths and travel times are incorporated into 
the computation of acoustic fields, the two approaches were 
seen to give very nearly the same result. 

Brokhovskikh [13] has shown that the lateral displa- 
cement A of a bounded beam for above situation is given 

A= - de/da 
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where €■ is the phase of the reflection coefficient given in 
(2.33) and is 

6- = 2 tan 

and a. is the horizontal component of the wave number and 

is given as oc — cos 0, where co is the operating 

C 1 

frequency of the source. 

The time displacement "t associated with beam displacement 
is given as 

-y __ d €: 

c “ ‘dm 


L 


p . , .2 

cos ©— ( c C 2 ' “{ 

" — ... — j 

1— COS 0 


Evaluating A and T , we will have 

2 _ 2 


2 ? 1 P 2 c 1 c 2 ^ C 2 “ C 1 ^ C0t 9 
A = ^cJ^^-T^P 7 sin 2 * 


-2 2 \ 

Pi c 2 5 


and 

r = 


Pi p 2 (c 2 2 “ c l^ 2 cote c ose ^ 

7(7 2 2 cos 2 G - c 1 ) 1 ^ (F/c/cos 2 © +P 2 2 c 1 2 sin 2 ©-F 1 2 c 1 2 ) 


: or evaluating A and T we make use of direction cosines 


given as 


w 


y = «*•-- sin© 
l c 1 


w / 2 0 \l/2 

Y 2 = c~ (cos y - if ) 


a = ~ - cos© 
C 1 


0 being the angle the ray makes with horizontal, 



e 


2 4 


Thus A is a function of frequency. The value of 
A is much less than a meter per reflection for frequencies 
that are of interest to , this thesis. For moderately long 
ranges, the additional path to be considered is negligibly* 
small and so the ray displacement contributions are ignored 
in further discussions. For example 
for f = 30 KHz 
we have A = 5 cm* 

2,4 MULT I PATH STRUCTURE OF THE CHANNEL: 

Up to this point in this section, we developed various 
equations to compute ray path parameters. In this section vie 
further investigate the nature of the ray paths. 

In Section 2.3 we have seen that for a medium in 
which velocity profile is a complex function of depth, every 
point on the ray has to satisfy the Snell* s law and a ray path 
is confined to a plane nor sal to xy plane. Thus the various 
ray paths can be traced in principle. Depending on the 
complexity of the profile and the permissible resolution 
between paths, for each ray starting from the source, vie 
consider the range over which the velocity 1 c* can be 
considered to be a constant or a linear function of depth. 
Accordingly we approximate the ray path either with a straight 
line or with an arc of circle. <.e trace the path considering 
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the successive segments over which velocity profile is 
approximated. Thus we can draw all the ray paths that 
eminate from the source. Depending on the location of target, 
we can compute the number of rays that intercept a target. 

For a given configuration of source and target locations and 
channel depth, wc may be able to plot in finitely many 
paths connecting the source and the target. 

For ease in handling, the following nomenclature is 
often adopted in classifying various rays. Among the group 
of rays that are obtained, we have a direct ray which is the 
line of sight ray. The remaining rays are ordered according 
to the number (n) of bottom reflections the ray has suffered. 

This number is called the order of the ray. Now n can have 
all integer values starting from zero. For n=0, i.e. for 
the zeroth order ray we have only one surface reflection. For 
each order of bottom reflection with n > 1, we can associate 
a set of four paths; one of them will be having n+1 surface 
reflections, another one will have n— 1 surface reflection and 
two will have n surface reflections. The direct ray, and the 
rays having order n=0,l,2 are illustrated in Fig. 2.6. 

we assume the velocity of sound c at a point in the 
medium is linearly proportional to the depth of the point in 
the medium. Depending on the positivity or negativity of the 
gradient, location of the source and the initial ray inclination. 
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each ray travels in one of the four characteristic ways. 

These categories axe: 

i) the ray which gets reflected at both the interfaces 

ii) the ray which gets refracted at both the interfaces, i.e. 
the ray will not reach the interfaces but curves away 
due to the linear dependance of c on depth. 

iii) the ray which gets reflected at surface and gets 
refracted at bottom (RSR) 

iv) the ray which gets refracted at surface and gets 
reflected at bottom (RBR). 

These four types of rays are illustrated in the 
Figure 2.7. 

M.J. Jacobson [14] considered the ray paths under 
RSR and R3i. categories and gave formulae to find the number 
of paths arriving at a given point from a source, path 
losses, path travel time and the initial inclination angles. 

We now consider the approximations we make and the 
simplifications that result on ray path tracing in the channel 
of our interest. 

We assume an iso-velocity shallow water channel 
which is having flat boundaries. Therefore ray paths 
travel in straight lines irrespective of the location of 
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the source and all rays belong to the category (i). The 
rays will have specular reflection at the boundaries. The 
direct ray will have spherical spreading and due to good 
reflectivity conditions assumed for boundaries, rest of the 
rays will have cylindrical spreading. Hence received energy 
decreases in proportion to the square of distance for 
targets at closer ranges and decreases in proportion to the 
distance for longer ranges. The under water weapons will be 
generally active in the ranges ranging from 500 meters to 
2500 meters. Therefore on the average we assume the loss is 
due to cylindrical spreading. 

From the discrete ray structure shown in Fig. 2.6, we 
can see that the direct ray arrives first at the target and 
is followed by either S ray ox n=l ray depending on the 
location of source and target. Then the higher order rays 
will arrive. As the order of the ray 1 n f increases, the ray 
inclination at the source will become steeper. If the 
initial inclination is steeper than the critical angle (© ) 
defined in sec. 2.3.2, theray suffers reflection loss in 
magnitude. The number of rays that are to be considered are 
decided by the receiver and transmitter beam patterns. Among 
them we may neglect the rays that have inclination angles 
steeper tfyan critical angla and undergo more bottom reflections. 
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On the average fox a typical values of c i» ^ 2 * c 2 a 

ray suffers approximately 10 dB energy loss per bottom 
reflection of inclination angle exceeds critical angle. 

Till now we consider the boundaries to be flat, but 
in a physical situation the bottom structure may have either 
positive slope or negative slope. The ray tracing can be 
carried out with a simple modification of the above procedure. 
This is discussed by G.V. ,-.nand and P. Balasubramanian in 
their report [l5j. 

2.5 EFFECT OF TARGET MOTION ON RECEIVED SIGNAL CHARACTERISTICS: 

Using ray theory, we have seen that shallow water 
channel gives rise to a discrete multipath structure to the 
received signal. The received signal energy now is obtained 
by superposition of the energies arriving along these paths 
after getting 'reflected* at the target. If 1 n’ significant 
paths are to be considered each of the path energies may be 
characterized by an attenuation coefficient (A n ), arrival 
time (T ) and initial path inclination angle (© n ) . 

We now consider the effect of target motion. A 
moving target introduces a doppler shift in the frequency of 
the received signal with respect to the transmitted signal 
frequency. If V is the velocity associated with the target. 
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the doppler shift is defined in terms of the component of 
the velocity along the direct ray. we assume direct ray 
path length is large compared to the target movement over 
the time interval of interest and hence the only effect 
on the received signal is changing of phase with respect 
to receiver time frame ' t 1 or equivalently doppler shift in 
the transmitted frequency. 


'lie may now define a doppler shift (k n ) along each 
path by considering the velocity component of the target 
along that path and it may be given as 


k 


n 


= 2 


5a tH 


where © n is the velocity component of target motion 
along nth ray. 

We now argue that difference between the doppler 
shifts of any two discrete paths is small as we view that 
the doppler shifts are due to the changes in the path 
lengths due to target motion. Hence v/e propose a common 
doppler shift for the total received signal which is inde- 
pendent of the order of the discrete ray path, we do this 
first by computing the phase of the received signal along 
each path and then showing that the change in phase of the 
combined signal is independent of the order of the constituent 
rays by using some valid approximations. 
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We first write an expression for the combined 
received signal in terms of it’s constituent paths and 
for a general target motion. 

Consider the geometry of the system as shown in 
the Fig, 2,8* 

Let the transreceiver (TR) be located at origin ’O' , 
Suppose that the target at any instant of time follows a 
path a with arbitrary velocity 1 c‘, sound velocity could 
be a function of time. Target depth may change with time. 
Let the instantaneous position of the target (T ) be given 
by P, Then OP is the slant range and is denoted as R(t). 
Consider the plane X*Z which contains the directed line OP. 
The projection of OP on OX 1 is OP’ and is called horizontal 
range. The target is assumed to follow a path o as shown 
and V(t) = is the velocity vector. We assumed that 
* n* discrete ray paths eminate from transmitter, get 
Reflected* at the target and arrive back at the receiver. 

We also assume that the return path from target traces 
the corresponding forward path. Let © n be the angle nth 
ray (PR^) which leaves the target, makes with OX' . The 
convention for the sign of angle we adopt isj © n is 
positive if the ray travels in positive direction of 
z-axis (i.e., towards increasing depth) anc! © n is negative if 
the ray travels towards decreasing z. 
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We have seen in Sec. 2.3.4 that the first ray to arrive 
at the receiver after reflection at target is direct ray. 

It is followed by rays having bottom reflections of order ’ n' 
zero and onwards. Since the target is moving, the instanta- 
neous position of P follows the path and accordingly the 
vertical planes containing different rays may be different. 

Since the target is moving, the signal associated 
with nth ray at the receiver time 1 t* corresponds to a pulse 
reflected at a target position at an earlier time. For nth 
ray, let the ray leave the transmitter at time * t 1 and reach 
the receiver at time 1 1 ! . Then we have 

t n = t - 2 T n (t) (2.34) 

where T n (t) is one-way travel time. 

Here we see t to be direct function of t and so we may 
write 

t n = F n (t) (2.35) 

Let the transmitted signal is rectangular pulse with 
carrier frequency 
write, 

s(t) = Pie [Exp (jut)] 0 < t < T (2.36) 

Now the received signal along nth ray is written as 
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* n (t)=Re[A n (t-T n (t))exp j(co Q (t-2 T n (t)+S n (t)) ] 

(2.37) 

where A^(.) is the attenuation coefficient along the path. 
This includes path losses and reflection coefficient of the 
target y n . S n (t) i s ^he phase introduced due to boundary 
reflections and also includes the phase {S n ) of complex 
reflection coefficientof the target. The nth ray may have 
N bottom reflections and m surface reflections along it* s 
path. Then 

$ n ( t) ~ N^( t) £ n ( t ) + 7t t) + (2.38) 

where- is the phase change per one bottom reflection* Each 
surface reflection introduces a phase change of n. In 
principle number of reflections and phase change due to 
bottom reflection could be functions of time. 

\(*) i s given as 

Afc(t) = jaj L^t) [B n (t)]^ n(t) (2.39) 

where I^^(t) is two way spreading loss and we assume the 
path length difference in forward and backward paths is 
negligible over the analysis interval. 

B (t) is the path loss due to one bottom 
reflection ja^j is the modulus value of complex reflection 
co-efficient. 
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Now we want to represent the received signal as a 
complex envelope with an associated carrier frequency a) Q . 
First we write 


r(t) = Re[£ A^t) exp j J o> Q t -2“ 0 T n (t)+ £ > n (t)}-] 


= Re[s A n (t) exp j (u Q t-0 n (t)] 

(2.40) 

where 0 n (t) = -2w 0 T n (t) + S n (t) 

(2.41) 

Then we write 


r(t) = A(t) exp j(w o t+0(t) 

(2.42) 

Therefore, 


A( t) = U Sin0 n (t)] 2 + [S A n (t) 

cos 0 n (t)] 2 


(2.43) 

and 


. 2 A^(t) sin 0 ( -b) 

0( t) - tan - 

2 A (t) cos 0 (t) 

(2.44) 


We now introduce some assumptions and approximations 
and derive the composite signal for the case of interest 
to us. 

2.5 ISO-SPEED CHANNEL: 

We assume that the velocity of propagation ’ c» is a 
constant throughout the medium. The channel is bounded by 
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horizontal boundaries. The target is assumed to follow a 
horizontal straight line path with constant speed v. 

As shown in Fig. 2.9, let P be the location of the 
transreceiver and Q be the location of target at t =0. how 
R Pw gives the slant range at t=Q. Let PC<* P^q de 
the horizontal range at t=G, and y be the angle which the 
target path projected on to the horizontal plane makes with 
PCJ at t=*0. IL is the instantaneous position of target and 
r!(t) is the instantaneous slant range. It’s projection on 
AY plane is R^(t). 

From Fig. 2.10 we observe, transreceiver 1 P’ is at 
a depth of H Q , target is at a depth of H T and the channel 
depth is H. Vertical height of PQ is ’n* . As we have 
seen previously, each order of reflection is associated with 
a set of four rays. Each ray is characterized by its initial 
inclination. If ©; n is the angle mode by the nth ray, we 
follow a convention and say 8. n > 0 and a — 4-1, if the ray 
leaves the source downward? else © n < 0 and a = -1, and 
also p = fl if the ray reaches target/receiver from above, 
else J3 = -1. This convention helps in calculating path 
lengths easily. Some typical ray oaths in a moving target 
environment are shown in Fig. 2.10. Here the rays may not 


lia in the same plane 
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First we calculate the ray path length (one way) for 
N bottom reflections assuming target is stationary. For 
this we can use the method of images of the source that 
contribute energy at the target. Each ray may be identified 
by considering the number of times it crosses the surface 
and bottom or their images. From the Fig. 2.11 we see that 

R no =» nth ray path length 

= [R 0 2 +(2m+aH R + j3 H t ) 2 j 1 / 2 (2.45) 

where N is the number of bottom reflections and ft is the 

o 

direct ray length. 

In Fig. 2.11, the line segment from pto 0 indicates 
the Oth order ray path (surface reflection only) and line 
segments from pto ^ denote the set of 1st 

order ray paths. Ray path lengths can be easily computed 
from the geometry. For example, Ray . corresponds a path 
having two surface reflections and one bottom reflection. Its 
length is given as £k o 2 +(2H + • 

Therefore equation (2.45) gives nth ray path length 

/ 

as F* for a target at rest. Then one way travel time 
is given as 



(2.46) 
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From Fig. 2.9 we can compute the path length fox nth ray 
for a moving target. 


= C R no + v \ 2 + 2vt „ «HO cos r ] 1/2 


= fc 2 t L + v2t n 2+2vt n s„U- rV ) 1/2 cosy ] 172 


no 

where h = H H. 

T 6 


n o’ 


K 


(2*47) 


From equation (2.34) we write now 


* n = t-2T„(t„) 


n rv 


= t-2c*' 1 [(v t n ) 2 ^2 R q v t n (l- ^-2 ) 1/2 

R o 


cosy 


+ c 2 T 2 jl/2 

no ■* 


— t—2[M 2 t n *'+2Fs. 0 c 1 M t n (l~ ^—^) cosy 

B o 


. t 2 -, 1/2 

no J 


(2.48) 


where M is called mach number and defined as 

M = ~~~ (2.49) 

Solving equation (2.48) as a quadratic equation we have 

2 

t n (t) = F n (t) = [t + 4 R o iw(l )^ 2 o' 1 cosy j 

R 
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- [4 T 2 ft + 8 R Mt (X- 
no ° 


R 


2 

2 


4M 2 {t 2 ~4T 2 +■ 4(R 2 -h 2 ) 

' no o ' 

x (l-4]\; 2 )~ 1 


J 1 / 2 c‘" i cosy + 
c ~ 2 cos ^)] 1 / 2 X 

(2.50) 


If we make realistic assumotion now like << 1 and H << R 

* i 

(i.e., velocity of target is much smaller than velocity of 
sound and depth of the channel is much smaller than range), 
then equation (2*50) simplifies to 

F n (t) = t + 4 R q M c" 1 cosy - 2(T 2 q + 2 K q h c“ 1 (cosy)t 

+ M 2 t 2 ) 1 / 2 (2.51) 

F (t) is the transmitted pulse time frame (*t 1 ) which 
n n 

has limits 0 < F n (t) < T where T is the pulsewidth. i;e 
expand t n (F n (t)) using kaclaurin 1 s series. 

, 2R k c"^cosy 

F n (t) = t-2T n0 + 4 R o a cosy 2_ t 

no 

(2.52) 

2R [A c~^cosy 

Consider t(l — — — - ) in equation (2.52), where 

1 no 

T no is the travel time of nth ray when target is stationary 

and it will be in the order of 2R Q c"^. Even for cosy = l,. 
2R M c’ 1 cosy 

is negligable compared to '1' if M << 1. 

no 
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Therefore, we write 


F n (t) * t~2T +■ 4R M c "* 1 cosy 
n no o 1 


(2.53) 


where t = o (2R /c). 

o 

Vie now proceed to find expressions for the attenuation 
coefficients, phases and inclination angles associated with 
various rays. 


The inclination angle 0^ for the nth ray is easily 


seen from the geometry of Fig. 2.11 as 
-I 


Q =* tan’ 
n 


[ 


P H t - a H r + 2 N h « 




J 


( 2 . 54 ) 


Spreading loss associated with nth ray is computed 
with the assumption that they on the average are cylindrical 
losses, we have for one way path from equation (2.28). 

L~^ 2 (t) = | nth ray path lengthj 


= c** 1 ! T 


2 

no 


+ M 


+ 


2R c*" 1 M(cosy) t I 


-1 /: 


n' 


(2.55) 

Let B n be the boundary reflection loss of nth ray 
per one bottom reflection derived from equation (2.31). 

Then for N bottom reflections along nth ray we have total 

n N n 

boundary loss = [B n J 
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Therefore total amplitude loss is given as 


N. 


S n (t) = [3 n ] " L- X (t) | a n [ 


(2.56) 


where is the complex reflection coefficient of the 

target. 

Now we linearize L^^t) as follows s 


-1 


L n * (t) = c~ 2 [l n l + ivi 2 F n 2 (t) + 2R 0 M F n (t) c" 1 cosy] 


-1 


-1 


0 Q M 2 F 2 (t)+2R M F (t) C -1 C0SY , 

= c" 2 T' 2 [1 + 2 - V ~ ' ] _1 

no j 2 

mo 

(2.57) 

Again using the procedure in simplifying (2.51) we have 


L n X ( = c "^ T no U" 2 ** cos ^ ( !" “ 2)) 


(2.58) 


where T „ is the order of R /c. 
no o 

Hence total loss A n (t) is given as 


N 


A n (t) = [B n ] n jaj (c T no )“ 2 (1-2 M cosy(|^ -2)) 


ct 

‘o 


(2.59) 


The total phase is given as 


S n (t) = [N + i-SjSi ] % + (2.60) 

where first term gives the phase change due to surface 
reflections and second term gives the phase change due to 
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bottom reflections and€. is the phase advance due to one 
botton reflection. Here we assume the number of reflections 
H and phase change 6 are not function of time as the order 
n of the ray is not large. Hence S (t) is not a function 
of time. 

Knowing the values of A^t), S n and T n (t), the 
received signal contribution from the nth ray can be 
computed. It is given as 

r n (t) = ReO n (t) exp j | o> 0 ( t-2T n ( t) ) + ] f 

= H«[A n (t) exp j £ w Q t + 0 n (t)^ ] (2.61) 

where, 

0 n (t) = -2“ 0 [Tn 0 +M 2 F n 2 (t)+2R 0 K F n (t)c“ i cos Y ] 1/ ' 2 

(2.62) 

Consider the expression given in square brackets in the 
above equation; and using the simplification procedure 
adopted in case of equation (2.51), we will have 

0 n (t) = -2w 0 [T n0 + M(cosy)t-2K 0 M c^cosy]*^ (2.63) 

How we write the total acoustic field received as 
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N. 


(t) = Re[ X [B n j' n jaj (CT no )~ 2 (1-2M cosy(|^ -2) 
n o 


ct 


x exp j 5 (o t - 2 w T -2to M( cosy ) t+4(«) R Mc~ i cosy+S ? ] 
L o o no o v ' o o ’ n- 

Rewriting the above equation as 


r(t) = Re$£[i -2M cosy(~i -2))exp j(-2w o M(cosy)t 


-1 


N 


+ 4 w 0 R 0 Mc '" C 0 $ Y)j x L £ CB 3 n |a n l(c T no ) 

n 

ju t , 

exp j(-2 Vno^n^ * e J j 


-2 


(2.64) 


Thus the received signal is seen to be the product of two 
terms, one of which is a function of the target velocity 
alone and the other - under the summation sign - is the 
composite signal corresponding to a target at rest. Denoting 
the later term by A Q exp(j0 Q ) we write 


N 


A 0 exp(j0 ) = E [B n ] ' n |a n | (CT r 2 exp j(-2w o T no +S n ) 

n 


(2.65) 


Now we have 


r(t) = Re| [( 1-2M cosy (g- -2)) A Q exp j(-2a) o M(cosy)t 




* J LU <* ** 

+ 4 w 0 R 0 M c" 1 cosy + 0 Q )] e ° j 


( 2 . 66 ) 



42 


Comparing with the above, the r(t) given in equation (2.42) 
we observe 

tf(t) = 0 Q - 2 w q M cosy (t-^c" 1 ) (2.67) 

and 

A(t) = A q 1-2M cosy (~^ - 2) (2.68) 

o 

•«e thus finally have, as M <<1, 

i j<*> t 

r(t) « Re [A exp i(0 -2w M cosy( t-2R c~ x ) ) ]e ° 
v o o o 

(2.69) 

or write 

x(t) = Re [A 0 exp j(+ upt)exp j(0 Q ~4w o M cosy R Q c’* i ) ] 

iw t 

e ° (2.70) 

where, 

Wp = -2 w o M( cosy)t (2.71) 

We now observe that the doppler frequency (up), we 
obtained in the analysis is same as the doppler frequency 
one obtains in classical sense for a situation where the 
target is moving radially. 

Rewriting equation (2.70) as 

r(t) = Re [A 0 exp j(0 o -4w o M cosy R^* 1 )] exp j(o) Q +w D )t 


(2.72) 
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In the above equation, the term, 

«*o - 4 “o M cos l r V ' 1 = - 2 “o T no _4 “o M cosy R o c_1+S n* 

Since M<<1 and T is in the order of R c”^, we 

no o 

simplify it as 0 o (O), 

Therefore, 

M 0 

i(t) = Re [S [B n ] )a n ] (c T n0 )- 2 exp(- 2Mo T no+ S n )J 

exp j(w Q + w^t (2.73) 

Therefore the composite received signal for a moving 
target situation is seen to have amplitude equal to the 
amplitude computed when the target is assumed to be at rest, 
and having a carrier frequency shifted by from trans- 
mitted frequency , where co^ is due to a moving target. 
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2*7* Results: 

We want to investigate the variations on time 
delay differences and attenuation coefficients j a^j 
in a shallow water channel with known source and target 
depths* 

The values of the parameters we consider below are 
typical to a situation where a torpedo is homing on to 
a submarine* The following parameters are assumed: 

/ 

i) Velocity of propagation in the channel : c^ = 1500 m/sec. 

ii) Velocity of propagation in channel sediment: C 2 = 1600 m/se*, 

3 

iii) Density of the channel : = 1000 gm/cm • 

3 

iv) Density of the channel sediment : - 1250 gm/cm • ' 

v) Channel depth : H = 300 m. 

vi) Source frequency : f = 30 KHz* 

vii) Transmitting pulse width : T = 100 msec, 

viii) Target velocity : V - 15 m/sec, 

ix) Direction of moving target : y = 75°, 

x) Target range : R * 500 — 1000 - 1500 — 2000m* 

xi) Target depth : Bp = 0 - 50 - 100 - 150 - 200 - 250-300®* 

xii) Source depth : % = 50 - 150 - 250. 

xiii) Beam pattern : Vertical : 40°} Horizontal : 80°. 
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xiv) Target strength : IF 70° < y < 90° : TS » 300 

else IF 0° < ft n < 10° : TS = 30 

10 ° < ®n < 20 ° : T S = 35 

20° < © n <45° : TS = 60 

45° < 9 n < 90° : TS = 150 

We compute the desired parameter values, for 
direct ray and the sets of reflected rays with order 
N — 0,1,2, 3, It is observed that for N _> 3, the inclination 
angle © n of the rays are out of beam pattern. 

The following parameters are computed for each ray 
and the equations which are used are shown in brackets, 

i) Critical angle (0 ) : from equation ^z.3a-> 

ii) Inclination angle (^) : from equation ( 1 -^) 

iii) Time of travel (T) : from equation 

iv) Spreading loss (SPLS) : from equation 

V) Boundary reflection loss(BLS): from equation 

vi) Bottom reflection phase <- : from equation 

vii) Total phase introduced due 

to boundary reflection : + 

where M is the number of surface reflections, 

viii) Array factor AF : using © n and Beam pattern specifications 

ix) Target strength: TS: Using y, © n and target strength 

specifications. 

x) Total attenuation coefficient j a\J : SPLS x BLS x 

AF x TS, 
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Then normalizing the parameters T and j ILj with 
respect to direct ray parameters, we compute time delay 
differences and normalized attenuation coefficients j fcuj 
as 



where and | 'a d | refer to the direct ray. 

The following results are observed: 

1. For a given range of the target and source depth, time 
delay difference (Z^) of any ray path is linearly 
dependent on target depth (since the linearized version 
of the equation (2.45 ) is used). The gradient of the 
ray depends on the type of the ray. Refer graphs 2.1 
and 2.2. 

2. For a given source and target depths, the time delay 
differences (Z^) depend non>»linearly on range of the 
target. Refer graph 2.3. 

3. For a given source and target depths; graph 2*4 shows 
the attenuation loss with respect to range. 

4. The spreads of time delay differences for various 

rays are shown in Table 2.1 The rays that have higher 
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order have larger spread; i.e., if Z ^ are arranged 
in increasing order, then we -observe 

± < £ Z. for *1. < Tj 

where & ^ is the spread of the delay 

This fact will jpe useful in Chapter 3, in computing 

the correlation coefficients and their probability 

densities, among signals which arrive along different 

paths. 

5, The number of significant paths that contribute to the 
composite received signal at each range are shown in 
Table 2*2# 

6# Assuming each range bin is of transmission pulsewidth 
’ T T , starting from direct ray time delay 2T d , the time 
delay differences are arranged in corresponding range 
bins. This observation is helpful to see how the return 
signal energies via different paths are distributed. 

Some typical results are tabulated in Table 2 #2. 

7* Table 2.3 gives the type of rays that contribute to 
received signal at each given range. We observe that 
the number of rays contributing to the channel are same 
irrespective of source and target depths. The only 
change we see for 1000 m and 1500 m., is due to the 
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the symmetry of ray paths with respective to source and 
target depths. This result vail be useful in estimating 
the arrival times as will be discussed in Chapter 3 , 

( 2 . 1 , 2 . 2 ) 

8, In Graph [_ we observe at a given range and source depth, 

the variations in time delay differences. The values 

taken for range (K)j source depth (H R ) and target depth 

(Hy) are arbitrary but typical, Vve observed that some 

time delay differences (t-2> ^5 > "*6 etc.) are increasing 

with respect to target depth whereas the other time delay 

differences etc.) are decreasing with increase 

in target depth. This implies that as H R changes some 

paths are moving away from direct ray and some are moving 

f or 

closer. Therefore £ a change in Hy, there is a net amount 
of change in this movement. This result will be of 
considerable importance to analyze the performance of 
conventional result in Chapter 3 . 
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Table 2.2: Total number of arrivals and their 

distribution with respect to direct ray 
arrival time for different ranges. 


Range 

(R) 

Vertical 

height 

(hi 

Total 

arrivals 

Arrivals with respect 
range bin 

1 

2 

3 

4 

5 

600 m 

< 50m 

2 

2 

— 

— 

— 

— 


>, 50m 

2 

2 

- 


- 

- 

1000 m 

< 50m 

4 

2 

2 

- 

- 

— 


> 50m 

4 

2 

2 

- 

- 

- 

1500 m 

< 50m 

6 

3 

2 

1 

- 

- 


> 50m 

6 

4 

1 


B 

- 

2000m 

0 m 

7 

m 

3 

H 

■ 

— 


50m 

8 

■ 

2 


H 

- 


100m 

8 

■ 

2 



- 


150m 

8 

■ 

1 

2 

i 

- 


200m 

8 

■ 

1 

3 

mmm 



(h) : Vertical height between source/receiver 
and target depths. 
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Table 2.3: The type of rays that contribute 
at each range . 


Flange 

U) 

Vertical 

height 

(h) 

hays that contribute 
at receiver 
(Ref. Fig. 2.6) 

500 m 

0 - 250 m 

1, 2, 3 

1000 m 

< 250 m 

1, 2 f 3 f 4 


250 m 

1, 2, 3, 5 

1500 m 

< 250 m 

1, 2, 3, 4, 5, 6 


250 m 

1» 2» 3, 4 f 5, 7 

2000 m 

f 

0 - 250 m 

1, 2, 3, 4, 5, 7, 9 


(h) : Vertical height between source/receiver 
and target depths. 


CHAPTER 3 


RECEIVERS FOR TARGET DETECTION 
3.1 PROBLEM FORMULATION: 

In Chapter 2, it was seen that the signal received 
from a shallow water target is a weighted sum of a finite 
number of replicas of the transmitted waveform with 
different delays and a common doppler shift. The number of 
paths that contribute significant energy at the receiver are 
limited by i) the beam patterns of the receiver and the 
transmitter and ii) the number of bottom reflections a ray 
.path undergoes before it reaches the receiver. i 

In this chapter different receiver structures are 
considered and their relative performances studied. All the 
receivers are aimed at the detection of a target under the 
usual binary hypotheses testing [l6]. I 

t 

Let the transmitted signal be a rectangular pulse with 
carrier frequency w o . Accordingly [ 

$ T (t) = f2 Re[y£ t f(t) exp( jo) Q t) ] (3.1) ; 

A | 

1 

where s 

f(t) = -Yj 0 ^ t < T \ 


0 


otherwise 


(3.2) 
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E^t is the transmitted signal energy and the complex envelope 
f(t) has unity energy. Following our discussion in Chapter 2, 
we write the received signal component arriving via the ith 
path as 

j ( CO -rUU. ) t 

Y i 5 (t) = Y2 &e [yu t | a.| f(t-2T i )exp j ( -^T . +S> . ) ]e 0 ^ 


(3.3) 


where, 

is the attenuation suffered by the ray along 
the ith path. 


is one way travel time, 

is the phase introduced due to reflections at 
boundaries and target, 
and Wq is the doppler frequency. 

The composite received signal may be written as; 


y( t) = f2 Re [VE t 


M , .. j(u n +OK)t 

E i } a j .|f(t-2T i )exp J (~ 2w 0 T i +S i) 3 e 


(3.4) 


where M is the number of significant paths. a^* s are 
modelled as complex Gaussian random variables, bach of 

constitutes two parts . ^ complex oe terminr s cic variaol 

a describes path losses and boundary losses. And a 2i . 


I 

t 


8 
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zero mean complete Gaussian random variable, describes 
the reflection coefficient of the target, which we assume 

to be a slowly fluctuating point target. a 0 . has variance 

2 ‘ _1 
2^-, and the phase is uniformly distributed over 

[0,27cj » 

Now we write 



The phase distribution of is seen to be uniformly 
distributed with a bias component which is equal to a^' s 
phase, iVe also note ja^l a function of target's aspect 
and altitude angles with respect to receiver; which in turn 
depend on path inclination angles Hence a^ have in 

general unequal variances. 

Ne now consider the time delay (T^) and attenuation 
coefficients (am) normalized with respect to direct ray's 
time delay (T d ) and attenuation coefficient (a^). That 


is define 
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ti < T i- V 

(3.6) 

r* 

a . 

b. A __JL_ 

1 ~ ^ 

(3.7) 


for i=l... M. 


T^'s are now time delay differences and b^ 1 s axe nor- 
malized attenuation coefficients. If direct ray subscript 

is taken as ’l 1 , then , = 0 and 2cr, 2 = 1. 

i Ojl 

Rewriting the composite received signal along with 
noise component which is considered to be a sample function 
from a zeromean complex white Gaussian random process, 
we have 

y( t) = f2 Re[J YB t E i b i f(t-t i )j exp^ tj 

+Y"2 Re [ w(t) exp( jw Q t) ] (3.8) 

As we observed in Section 2.6, the composite return 
signal is having a common doppler shift. Therefore with 
no loss of generality we equate to zero. 

The complex envelope of the received signal is 
obtained under binary hypotheses testing as, 

y( t) = YE t b A fCt-rp + w(t) : H x , -,) 

w(t) : H q 



where 


0<t<T* 11 T + Max( j) ; j [1...M] 

and B[w(t) 'w(u)] = 2a 2 £ (t - u) 

hr 

Vi/e now consider various receiver configurations and 
attempt to evaluate their performances. 

3.2 OPTIMUM RECEIVER : 

The complex envelope of the received signal is given 
by equation (3.9) and is reproduced as 

y(t) = y£ t b i fCt-fj.) + w(t)i h x 

w(t): H 0 

where attenuation coefficients ( b^) are modelled as complex 
Gaussian random variables. Time delay differences ('tb) are 
functions of target range, target depth and source depth. 
Hence s can be seen to be some functions of target depth 
if range and source depth are known. Since target may lie at 
any depth in the channel, may be modelled as random 
variables, in a range ^ ^ A ^ *" L imin * 

Now we are interested in designing an optimum 
processor with above specifications, i.e., the problem is 
to observe y(t) over an interval [o,T‘ ] and decide whether 
H 0 or is true. The criterion may be either 3ayes or 
Neyman-Pearson criteria. 
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Our observation is a time continuous random waveform. 
The first step is to reduce it to a set of random variables. 

i.'e observe in equation (3.9) that the received signal 
y(t) is a zeromean random process, with variance K^(t,u) 
given by 

iq(t,u) = E I 20. (u-tp]+ 2o^(t-u) 

r i 

(3.10) 

'It e want to expand y(t) by using K.L* expansion method to 
obtain a set of uncorrelated random variables. 


We write 


Y ( t ) = 

where y^ = 

and £ 0 i ( t)5 is 
Expanding y(t) 


n 0j.(t) (3.u) 

/ T Y(t) 0*i (t) dt 

) 

a set of complete orthonormal functions, 
by using the above functions we will have 


Vi = W VE t P ij + •• H 1 

T 1 - • 

where p ^ ^ == / f(t-T'j) 0^ (t)dt. 

w. = / w(t) 0^(t) dt 

1 o 


(3.12) 

(3.13) 


and 


(3.14) 
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0j_( ! i = 1»2 ,..,n, has to satisfy the following integral 

equation: 

T 

^i ^i^) = / Kv (t,u) 0.(u) du 

o r 

where X ^ are eigen values* 

The above integral equation may be written as$ 

T * 

**!<*) = E t E l 2 a b E t / f (u-r.)]] 

J j O J J 

X 0 i (u) du + 2 crj 2 * 

Instead of solving the above difficult integral equation 
with a Kernel involving the statistics of the various delays, 
we explore the following alternative approach. 

Towards that end consider the expansion of y(t) in 
terms of complete orthonormal set o-f- functions ( t)3 , the 

«C. 1 

first M of which are given by f ( t-Th ) for some suitably 
chosen set of values £. ; i = 1...M, 

We define, now, two vector functions as follows: 

F(t)4 [ftt-tp f(t -t 2 ) ... f(t-^)] T (3.15) 

F(t)A [fft-q) f(t-r 2 ) ... f(t-r.)] T (3.i6) 

Then we expand y( t) over F(t) and the functionsjpgtjJ” . 

•.Ve obtain a set of random variables y^. 
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T* 

Yi =* S y(t) f*(t-r.)dt 

o 1 

= 2-j VE t bj p .(i) t.j.Bj 


"i : H o 

(3.17) 

T 1 

where P ( i) = / f(t-t.) f*(t-£.)dt 

J 0 j 1 

(3.18) 

T f 

and w. = / w(t) f*(t-E)dt 

o 1 

(3.19) 

For i = M+l,...; y^ are same for both the hypotheses 

and provides no additional information and hence can be 

ignored. Observe that vin are uncorrelated random variables. 
The cross covariance of y i and y^ are given as 

’^CVi Y k ] = E C r j(i) •' j( K )J + 2o w 2 iK 

(3.20) 

Nov/ we have a set of random variables y^, which are 
functions parameters b. and (°.(i). he want to now formulate 
the problem as a composite hypotheses problem, he first 

introduce some notation. Let 

ft £ tYp y 2 ••• V M ] 

(3.21) 

T 

B-^ Tbj_» ^2 *** 

(3.22) 

w 4 [w x , w 2 ... w m ] T 

(3.23) 
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Accordingly we will have 
B[b] =4 o 

and E[ B B + ] 'b (3.24) 

where A ^ is a diagonal matrix with diagonal elenents 

equal to 2a. 2 . 

D i 

We also define the correlation matrix, 

T ' 

/ F(t) F + (t) dt (3.25) 

o 

X * 

where P.(i) = / f(t-T.) 'f*(t-T / .)dt 
J o ^ 1 

The composite hypotheses problem is now to compute 


A(r ) = 


B. p P CR/B, P, H^pC'e ) p( J> ) d Bd ‘ 


P (R /H 0 ) 


( 3 . 26 ) 


Towards evaluating A (R ), we first consider evaluating 

A( r/p) by averaging A (H /S , P) over the space of B 

parameters . A (R /P) is a complex Gaussian vector under both 

H. and H, and hence its solution is obtained by computing 
0 1 

the conditional likelihood ratio, which is given in 


general as 
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A ( R/P) = 


(T ^) 1 / 2 exp[- | ( R + - M x + ) Q x ( R - M x ) 
O ^) 1 ^ 2 exp [- |( R + - M q ) Q q (R - M Q ) ] 


where, 

H[ R/Hj] = M x ; H[ R/H 0 ] = M 0 
£[( R-M 1 )(R + - M 1 + )/H 1 ] = K x 


(3.27) 


E[(R - M 0 ) (R + - M 0 + )/H 0 ] = K 0 

' y l - V 1 » Q o = V** 

In the present problem, we have 


R = VE t P T B + 7/ : H x 


S3 

W : 

H 0 

(3.28) 

a nd 

o 

il 

o 

IS 

II 



K i 

= E t p’b 

P + 

(3.29) 


= K-w = 2 % 2 

I 

(3.30) 


Then the conditional likelihood ratio is simplified to 


/\ ( k/P) = | r^| l/2 ex Pt| R + (Q 0 - } R 3 

K 1 


(3.31) 
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Now we need to compute to obtain a simple 

structure for the above ratio. 

Consider the signal component of covariance matrix 
K^» and it is given as 

I<is£. £ t pT A' b P (3.32) 

This is a non-diagonal matrix implying that the compo- 
nents of vector It are hot independent, iVe now diagonalize 
the matrix k^ by solving its characteristic equation and 
obtain the eigen values and eigen vectors. If K. ^ is not 
positive definite, we augment the eigen vectors with some 
more to form a complete set. We have seen K w to be a 
diagonal matrix with equal diagonal elements (2a^ ) and 
so the diagonalizing matrix (0) will not effect the matrix 
K ,-j. So write 

k- s = 0K ls r 

and + k^ (3,33) 

If 2cr 2 are the diagonal elements or eigen values of K.^, 
they represent the energy components of received signal 
along different basis functions and also their sum is 
always constant irrespective of the choice of 0. Now 

' - Q * is obtained as 
X o 
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Substituting the above matrix in conditional likelihood 
ratio, we obtain, 


(3.34) 


-r— 4 

A( R'/P) = ] 1/2 exp -r_ (3.35) 

1 °i +0 «> 2 ° u 


E i 


a. , „ 

_Jk v 2 

2 . _ 2 Y i 


a. 

1 0) 


where R* = i# 4- R 


(3.36) 


The above ratio is developed with apxiori knowledge 
of P i.e., P.(i) for i ± ^ - 1...M. Now if we have the 
probabilistic knowledge of P, then the optimum processor 
would be 

(ft* ) = / ( r'/p) p(p) dp ( i/ 37) 

p 

This integration in general will be difficult and we 
consider a case where M = 2. 
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where K uS ^ A A E t ( <^2 p 2 ^) + p|( i)) 

KazsA 6 ^ E t‘°b2p^ (2) + « b ^2 (2) * 

and ^is = t 'i2&^ C — S t^ a b 2 ]/ 2 ^ + °b 2 <°2^ 2 ^ 

Solving the characteristic equation we will have 

2 c 2 
^ 1 * A + J~B and " 2 = B * A^B 

if c 2 < (A-B) 2 /4. 
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then the conditional likelihood ratio is given as 


A( R * /p) = 



1 2 2 V C i + V 





and then 


/*(&') = / A( R*/P) P(P) dP, (3.37a) 

where p. (i) = j f(t-7“. ) f*(t-Tjdt. 

■J o 1 1 

Now p.( i) is a random variable and it's probability 

•J 

density depends on the density of Tr.. From ■ .(i) we now 

w/ 

have to compute joint densities 7^ an d /)2 to evaluate the 
integral (3.37^. 


This approach of evaluating likelihood ratio involves 
extensive computation and the receiver implementation would 
also be difficult} even for the case M = 2. 

But however, if we could compute the joint densities 
of and ^ and ass,ame high signal to noise ratios, we may 

be able to obtain a simpler structure. The eigen values 
and y\ 9 are the received signal energy components and if 
> i >> a w 2 , in equation (3.35), we see that 


A ( R* ) = 2 a o 2 exp [-^ J ( V ' 2 + Y ’ 2 )] £ *) ( 3 . 38 ) 

w 2 a 2 

w 
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where 


a 



P(/\ 1 ^ 2 ^ d ^l d ^2 


( 3 . 39 ) 


Likelihood ratio is obtained as. 





* 


The performance of this receiver can be easily obtained as 
t (R* ) is having central chi-square distribution. 

V;e can also observe that for high order M, the 
sufficient statistic is still same i.e., 

h r') = s “ =1 Y ’ 2 >< Y* 

i 

if high SNR condition is satisfied. But we need to compute 
Mth order joint density function and then evaluate the 
corresponding integral. 

For low signal to noise ratios, the receiver may not 
be as simple as shown above. By appropriate modelling of 
densities ofTp we may still be able to realize the optimum 
processor. 

\jje now consider some other receivers wherein we have 
precise information of the time delay parameters. 



77 


3.3 OPTIMUM RECEIVER WHEN TIME DELAYS ARE PERFECTLY KNOWN: 

In Section 3*2 we have seen an optimum receiver 
for nutation when B and p are random vectors and matrices 
respectively. This was found to be a difficult problem even 
for M=2, Now we consider an optimum processor for the 
situation in which time delays^. and therefore p are 
perfectly known. The performance of this receiver if 
evaluated gives a bound on the performance of the receiver 
discussed in Section 3,2. In Section 3.4 when we discuss 
a conventional receiver with perfectly known delays, the 
performance of the optimum receiver will give us an estimate 
of the degradation that will occur with a conventional 
receiver. 


We as seen in Section 3.2, again consider a set of 
orthonormal functions^ t) 'j wf which the first M functions 
will be F(t) (eqn. (3.15). This choice of F(t) is essential j 
to obtain uncorrelated noise components when y(t) is expanded.’ 
Under known P, y^ are zeromean Gaussian under both the j 

hypotheses, but will be non-diagonal. Solution of it’s j 
characteristic equation gives a transformation matrix which i 
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This is to be compared against a threshold y* where 

2 2 

Y* = jFlnr] + In — ' ~ ^ ~ (3 

CT vj 

to 

The receiver may be implemented as shown in Fig. 3.1. 
3.3.1 Performance Evaluation: 


In the hypotheses testing problems we know that the 
log-likelihood ratio is given by ( R) and is compared 
against a threshold for detection of target. i{ R) is a 
random variable and it*s density depends on the hypotheses 
considered. If the densities of £ ( R) are known under H 0 
and in a binary hypotheses problem; we can compute the 
false alarm and probability of miss obtained for a given 
threshold and signal to noise energy. 

As seen in eqn. (3^ro), the random variable t ( R) is a 
sum of squares M Gaussian random variables with unequal 
weights. 


The computation of py H ( L/H^ is a difficult job 
as it involves unequal weights for the variables. Therefore 
we compute the performance bounds for the above receiver 
through moment generating functions [l7j. These bounds 
are called chernoff’s bounds. 


These bounds are found by computing a function 
M- ( s ) given as, 
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^i(s) 


oo 

ln 1 lp y/\ (R/H 1^ S [p T/H 0 ( ' fi/H o)^' SdEl 


(3.42) 


where r s’ is taken as a parameter such that s > o. Then 
false alarm is given as [i7], 

Jl 


Pp S. Exp[(p-/s)-s p,(s) + %- p.(s)) erfc*(sVp,(s)))lj;' 


*S > o 


(3.43) 


and probability of miss is given by 


2 

P M exp( p.(s)-( l~s) p(s) + ^(s)) erfc*((l-s)f ia(s)) 

s < 1 (3.44) 

where p.(s) and p(s) are the first and second derivatives of 


p(s) 


The minimum bound is obtained by computing s min > 


given as. 


M-C s) i 


= T) 


s=s 


min 


(3.45) 


where r) is the threshold. For equal probability of occurrance 
and equal cost functions, we have T) = 0. 

Computing a satisfying the above equation, now 

false alarm and probability of miss are given as$ 
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,2., 


P F < Exp [p(s) + p(s)]eri c*(sfpXs ) ) | ."*.4.*,) 


S~$ . I- 

mm 


2 

P M - t^( s ) *£( 5 ) J erfc*( sfp,'(s) ) | ^..47) 

where 0 < s < 1. 


s=s . 
ram 


If s min^( s min ) > 3 » then 


erfc* (s f p(s )) 

mm mm fu(s . ) 

y mm' 


Similarly if (1-s . ) Vu(s . ) > 3 then, 
7 x mm Jr ' mm ’ 


~s 2 . / 2 

> mm' -•/ n 
^ min' 


erfc* ( (1-s )fp(s)) j 


exp ( li= s L .'£(s)) 


s=s min ^*(1 -s)VV(») s=s m in 


Therefore, 

p < i exp [p(s)3 j 

f2ns^ ’il(s) s ~ s min 

> 

p < i~. exp [p(s)] | (3.49) 

V"2ix(l-s) 2 p(s) min 

3 , 3.2 To Compute Performance Bounds for Equation ( 3 . 4 -o ): 

’We have seen in Section 3, 3.1 that p(s) is given as 

CO 

H(S) -In / [A(R’)] S P y ,/h ( R/H 0 )<JR'. 

^00 1 ' 0 
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We know A (R«) and p y ,/H o ( R/H c ) where A (R* ) 
.s given by eqn* (3«3 C ! ) where — q * is given as 

7~ 9 

i 0 1 

0 


r < - 4 - 




cr. 


JL 


,2 2 X 2 
{ a . +cr ) cr 
v M w ' w 


Then p(s) can be evaluated and is seen equal to 


M 


p(s) = in / [ % ( 


2 — 2 )] 

- 03 i=l a, +0 


Vn ^/2 M 


-s 


1 0 ) 


_ 71 [exp [| + 2 

i=l a. 


w 


vi 2 


,2 


■ TTT2 ^ X ,7T"~2 exp[ ~ 772 J dy l (3.50) 


a i + % 


]f2na 


n 


2cr, 


0.) 


Completing the squares in exponent and solving the integral 
p(s) is simplified as, 

, ( d 2 ) 3 ( c , 2 ^, 2 ) 1 - 8 

lit.) - I 4 in t i J 

d- s )°i + 


Now to compute Chernoff's bounds, we proceed as 
follows: 

i) Equation (3*S© ) is computed for various values of ' s‘ 
in the interval [0,lj. 
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ii) The minimum value of ja( s ) is found and corresponding 
’s' is the s min and at this 's' we have p(s) = 0. 

iii) We now compute p(s) as 


p(s) 


X , r ( °i /V> 

2 i L 77 r,v 2 


(l+(l-s) c.7a; 


iv) Calculate u( s . ), 

min' 

v) Define x A s min fV(s min ) and y* U-s min )lT V(s min ) 
now if erfc* (x) > 3, the false alarm is given by 
eqn. (3.48 ) or else by eqn. (3.^6> ). Similarly if 
erfc*(y) > 3, the probability of miss is given by 
eqn* (3.4x-"i ) or else by equation (3. *. 7 ). 

vi) The total probability of error is given by, under equal 
likely conditions, P r (e) = 2 P F ’ + 9 P IviA* 


The results on these bounds are discussed in Section 3.5. 
We now proceed to consider a much simpler processer called 
conventional receiver. 


3.3.3 Selection of Transmitted Waveform to -Reduce Receiver 
Complexity: 

The complexity of the receiver structure considered 
in Section 3.3 arises from the fact that the signal contri- 
butions reaching the receiver via the various ray paths, are 
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not orthogonal. One way of making them orthogonal would be 
to choose a transmitted wave form f(t), such that its 
correlation function 

T 

R f (r)A / f(t) f"(t+r)dt (3.51) 

o 

satisfies the property, 

R f (r) = 0 for all f > T c 

where, 

T < (3.52) 

ip 3 

It may be recalled that all i^j , are non-negative 

differences in delays, between the reflected ray paths 
and the direct ray. 

A linear frequency modulated waveform or a Barker 
coded sequence of pulses or a pR3S sequence, with appropriate 
choice of parameters will approximately satisfy the condition 
given in equation (3.52 ). 

Assuming as in the previous section that the time delays 
T'- are perfectly known at the receiver, with a transmitted 
waveform f(t) chosen as above, we can easily show the simpli- 
fication that will occur in the receiver structure, compared 
to the receiver shown in Section 3,3. For this purpose 
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instead of considering the subset of the orthonormal functions 

» i=l. . .M as the first M members of the COM set 

' "V i ' b . as tn Section 3.2, we now choose them to be 
' •'1=1 

> i=l...M. With such a choice in view of equation 
(3.52 )» it is easily seen that 

pjU) = 1 if 

0 otherwise (3.53) 

As before the hypotheses testing problem simplifies to 


Yi = 5 VE t bj PjU) + w i S H x 


w. 

1 


H, 


(3.54) 


The covariance matrix of the measurements is now a diagonal 
matrix with i,jth element given as 


Hy 


i T i 


] = 


2Ej 


°b + 2a < \ 
D. CO 

1 



Proceeding as in Section 3.2, the likelihood ratio is given 
as 


m 

A ( R) = n: [ 


M c 

—5“ p J exp C 2. “2 

i=l E + < +o,A 1 2a;(E t a"+a 


E a 2 

t a bi 2n 

2772 • N 


'tHoWWo 


to 


The log-likelihood ratio is; 


(3.55) 


£( R) 


E t ab : 


1 ( E t< + 


Yi l 


(3.56) 
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and the threshold y 


[inn - In 


M 

% 

i=l 



E t-V 


to 



The receiver implementation is shown in Fig. 3.2. 


The performance bounds of this receiver can be computed, 
his can be easily carried out with a substitution of 


Li .a. 
t b 


/( E t a b 2 + ° 0J 2 ) in Place of cm 2 / ( cr^+cr^ 2 ) in eqn. (3 . 50 ) 


t b. ‘ “00 
Therefore p(s) is given as 


< E t°b 2 + > 


2 ,1-s 


p(s) 


l In 


( 1-s ) £, a. 2 + a 2 
v ' t b^ a) 


and 


a* /a, 


p( s ) 


. 't u b. ' -co 0 

1 2 [ i nT—o ] 2 

l+(l-s)(E t a b 2 /a/ 


We can compute Chernoff bounds by computing s min 
and then substituting this in corresponding equations. 

The results obtained by this method are discussed 
in Section 3.5. 

3.4 CONVENTIONAL RECEIVER: 

Consider the equation (3.9) in which we assumed b i are 

zero mean complex Gaussian random variables . s, the time 

delay differences are assumed to be known. Depending on the 

values of T. and^i we have and overlapping. 

1 J J 
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The conventional receiver is designed with the 
/iew that only career ray and the white noise component 
are present. ihen knowing the direct ray time delay (T ), 
we compute 


Yi = / y( fc ) f (t - T. ) .dt 

n ± 


(3.57) 


or equivalently when we assume = 0, we have 



y(t) 


* 

f (t)dt 


Si Y2 t b. /0.(i)+ on s Hjl 


w i : H o 


(3.58) 


where, 

i(l) = J f (*" i> f *< t " i) dt (3.59) 

and-.(i) -are known. 

Therefore y^ is a sum of utmost independent zero mean 
complex Gaussian random variables along with a white noise 
sample, under Hence y^ is a complex Gaussian random 

variable under with zero mean and variance ) where 

>! 2 ■ ti E t 2 i (1) (3 ‘ 60) 

2 

Under H„, y, is having variance 2 a . Then 
OX ^ 

likelihood ratio is; 
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A ( 'R x ) = — ■ d - : 2 " exp t(Yi)' 


a l + a co 


"o' 2 *2 -3 £ 

a AaA+a z * 

CO v 1 (.0 


and sufficient statistics is 


t( R) = IyxI 2 ? y * 


where, 


y* = [lnr] - In 


co 

‘ 2 S' 

°1 +c u, 


„ 2 , „ 2 2, 
0 ( ff, +0 ) 

CO v 1 CO ' 


Performance of conventional receiver ; 


We have 


l ( R) = |h.. I 2 = |/ y(t) f*(t- t )dt| 2 ( 


and we know the receiver performance is given by, 


A 


EURp 2 /^ E [( R 1 ) 2 /H Cj ] 

E[| Rj 2 /^ ] 

2 E t S 3 . i (1) 

r 3- from equation (3 .60) 


0) 


Mow vue know that false alarm Pp and probability of 
Pq are related as; 

P F = (P D ) 1+A 


3.61) 


3.62) 


3.63) 


3.64) 


3.65) 


detection 


(3.66) 
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Now for various values of SNR., the performance of conven- 
tional receiver can be computed. In Section 3*5 we compare 
the performance of this receiver with the performance bounds 
obtained for the optimum receiver with known time delay 
differences considered in Section 3.3. 

3.5 PERFORMANCE. CALCULATIONS OF CONVENTIONAL RECEIVER; 

Th© performance of a conventional receiver can be 
evaluated from the plots of probability of detection (P Q ) 
with respect to SNR for given values of probability of false 
alarm (Pp). The relation between Pp and P^ is given as 
(from equation (3.66); 

Pp = (P D ) 1+ A (3.66) 


where , 




P± 2 u) 


2 

Here a. is the variance of the reflection coefficient which 

i 

also includes path and boundary losses. P^(l) is given as 

PAD = / f*(t-r 1 )dt 

x o 

p 

He observe that 2E^ cr^ to be the total energy that arrives 
via ith path and therefore 2E^ o b 2 p i 2 (l) is the amount of 
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energy that ith path contributes in the performance factor 

1 A’* is the signal to noise ratio, then 

2 2 

2 4 cr^ 7(1) gives the multiplicating factor, which is the 
i 

improvement in the amount of energy 

that arrives in range bin 1. This multiplicating factor, we 
call as multipath effect (M.E.)* Therefore, 


l\ \ m iZ » 


°b 2 Pi 2 ^) 


l 

As we observe from equation (3.66 ), the magnitude of M.E. 
.has direct bearing on HOC curves, computed for a receiver 
which receives only direct path. As M.E. > 1 always in a 
multipath environment, the receiver perforins better than the 
receiver which is operating in single path environment. 

Hence in the following we explore the characteristics of 
channel which effect M.E. and the performance of the receiver 




Define ’ h' to be the vertical cistance between source ; 

i 

and target i.e., h & where H R is source depth and Hj I 

is the target depth. j 

1 

The magnitude of M.E. varies with respect to range, j 
source and target depths, we plotted various graphs 

| 

depicting the variations of M.E. with the above parameters. i 
I've infer the following facts. I 
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~E?J£5LiL.i : This graph is plotted for the variations in M.E. 

with respect to range for different h, with H p = 150 m. If 
the source is located at the middle of the channel, the M.E. 
linearly increases with range irrespective of the target 
depth. But for a given range, the M.E. increases with h. At 
h = 0, the M.E. is minimum. This is due to less number of 
arrivals' into first range bin (RBI) compared to h^O (refer 
Table 3.1 ). 

he observe that the number of arrivals in RBI decides 
the initial value of M.E. and M.E. increases with range and 
this is because the net amount of overlap in RBI by other 
rays, is increasing. 

Gr aphs 3 .2 to 3.3 : This graph indicates the changes in M.E. 
with range, considering ’h 1 as a parameter, with H p = 50, 25Qm. 
If h is a maximum, i.e., in a situation where source is close 
to one boundary and the target is close to another boundary, 1 

i 

the M.E. is negligible at short ranges and has steep gradient | 

f 

with respect to range, he also observe, with a decrease in I 

h the gradient also is reduced. For h equal to 200m, the 1 

i 

increase in M.E. is 20 OA over a range of 1500m, and for j 

h equal to zero, the increase is just 10/- for the corresponding; 

i 

I? 

range. But it is seen that, for h=0 at short ranges, the [ 

M.E. is 70/» higher than the M.E. at h=200m. These facts can be ! 
explained as follows. 
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For h = 200m, in Table 3.1 we observe that three 
new paths arrive in RBI when target range changes from 500 
to 1000m. This contributes to the steep gradient corres- 
ponding to , this curve. Similar effect, but to a lesser 
extent is seen for other f h 1 and hence their gradients are 
less. For h=200m; there are no new arrivals after 1000m, 
but because the net amount of overlap increases in RBI with 
respect to range, the M.E. also increases. For h = 0 m., the 
net amount of overlap is not changing markedly and hence 
there is no relative improvement beyond 100m range. At 
range of 500m, the M.E. increases inversely proportional to h, 
because the amount of overlap of the two rays present in RBI 
increases with respect to h. (Refer Section 2.6). 

Gra ph 3,4, 3.5 : This graph indicates the variations in Ivi.R. 
with respect to target depth for a given Hp^, taking range 
as varying parameters for different curves. 

When the source is at the middle of the channel, then 
M.E. is almost constant with target depth irrespective 
range. These curves only show a small reduction in M.E. when 
h=0. This is due to number of arrivals at h=0 are less than 
number of arrivals f or h > 0 (refer Table 3.1 ). If the 
source is not at the middle of the channel, M.E. increases 
with h. For range 1000, the k.E. initially decreases with 
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h f because there are no new arrivals and the only one ^ 
is present and this increases with h (refer Table 3.1 ). For 
range 500m, the jvi.h. decreases inversely proportional to h, 
because there are only two rays in RBI and as 1 h' increases 
the second ray moves out of RBI. 

Now we propose the following conclusions; 

M.E. is a fs notion of both target range and *h' ; in 
other words M.E. is a function of number arrivals in RBI and 
the amount of net overlap in RBI. he also observed from 
equation ( 3.66) that P^ is directly proportional to M.E., for 
a given and SNR. We now propose that if the source is in 
the middle of the channel, the M.E, is minimum and therefore 
the improvement on P^ is also proportionally less. If the 
source is not in the middle of the channel, the detection 
probability is markedly larger at long ranges as h increases 
and reverse the case for short ranges. Fox a given source 
power, SNR decreases with range and we can effect a reduction 
in source power to an order of 6 dB if the target is detected 
with source located at maximum possible h at long ranges. 

For short ranges irrespective of range, we have approxi- 
mately 2 dB reduction due to M.E. 
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3.5.1 Comparison of Conventional receiver performance with 
performance bounds of optimum receivers in Sec. 

3.3 and 3.3.3: 

We computed the performance bounds of tho receivers 
suggested in Section 3.3 and 3.3.3 i.e, for receivers where 
delay times are perfectly known. The receiver in Sec. 3.3 
is having rectangular transmitting pulse, where as the one 
in Sec. 3.3.1 is having modulated transmitting pulse which 
will result in mutually orthogonal signal waveforms along 
different paths. The performance bounds are calculated 
following the development in Sec. 3.3.2. For the same false 
alarm the probability of miss for conventional receiver in 
Sec. 3.4 is also computed. The results are tabulated in 
Table 3.2. The parameter values used are arbitrary but similar 
results hold good for all sets of parameters of interest. 

The following inferences are made: 

1. The degradation in conventional receiver cannot be exactly 
found as the performance results corresponding to the 
optimum receiver are bounds. For a given Pp and SNR, the 
of conventional receiver is always less than that of 
optimum receivers. Hence the conventional receiver in 
Sec. 3.4 offers satisfactory performance. 

2. The conventional receiver designed for the receiver 

discussed in Sec. 3.3.3 will offer no improvement over the 
conventional receiver designed then the channel is thought 
to be unbounded. 




* 3k« tob<4. U bracket* indicate 
ttte Vlutabe*. 6$ QMlmU 
iS0U4ce U at: ft* wL<jdile o ft* Cfiavmel 





HulU^Afi 






H «.!&*><* IT €i$edtr 




0 SO 

190 

ISO 

Xoo«rs^- H ca 5o 

^"■'y 

'■2y.&' : y 

50 


. . f 

Sva^K 3 >4 








CHAPTER 4 


CONCLUSIONS 


1 In this thesis we made an attempt to model the shallow 
water targets, and then we investigated some receiver 
s oructures for detecting them under white noise conditions 
by an active SODAR. 

Towards this, we have used ray theory to develope 
the acoustic transmission in shallow waters and it resulted 
in a discrete multipath structure. Each ray path is 
characterized by parameters; travel time, spreading loss, 
and boundary loss and phase shift. Hence the received 
signal energy is the superposition of signal energies 
arriving along different paths after getting reflected from 
the target. He have seen the number of ray paths that 
contrioute significant energy at the receiver are limited by 
SODAR beam pattern and critical angle associated with bottom 
reflection. We then analysed the variations in total loss 
suffered by a ray and it* s travel time with respective 
different ranges of target from receiver, for source depth 
and target depth. 

For an iso-speed channel, assuming target velocity 
to be much smaller than sound velocity and channel depth 
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to be much smaller than direct ray range, we obtained the 
linearized model for the received signal for a SODAR from 
a moving target. The received signal energy was interpre- 
ted to be same as that of a stationary target with doppler 
shift in carrier frequency 'independent of number of rays 
that contribute at receiver. This common doppler shift 
for the composite signal allows us to look for a simple 
receiver configurations. 

The detection problem for the shallow water target by 
bODAR is formulated, be assumed the combined effect of 
path losses, boundary losses and losses due to reflection at 
target as attenuation coefficient of the signal arriving 
along the path and it is modelled as zero mean complex 
Gaussian random variable. If the arrival times are also 
treated as random variables, then the optimum receiver for 
this characterization is practical only at high signal to 
noise ratios. Assuming perfect knowledge of the time of 
arrivals, the optimum receiver was investigated and analysed. 
The conventional receiver for this characterization found 
to, be adequate. Usage of special signal waveforms (Barker 
coded, LFivi etc.) for transmitter pulse resulted in simpler 
optimum receiver structure. But the conventional receiver for 
this specifications gives no additional advantages. 
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Suggestions for Future Work: 

As an extension of the work done in this thesis an 
appropriate model for the channel where the sound velocity is 
a function of depth and the boundaries are of complex 
structure. The target may also move along any specified path. 

The model may also include the effects of reverioaration 
which then will be more realistic. An appropriate model 
for the arrival to times may be attempted to obtain an 
optimum receiver in a composite hypotheses situation. 

For the optimum receiver analysed in Section 3.2, we 
assumed to have perfect knowledge of time delay differences. 
Some efficient estimation methods may be investigated. One 
possible way of estimating travel times is as so follows: 

With an apriori knowledge of approximate range and 

receiver depth, the time delay differences (ts) may be 

modelled as uniformly distributed random variables. If 

A .= T . If T is the correlation time of the 

> r max i mm. o 

transmitting signal and KT Q = , then using a scheme shown in 

Fig. 4.1 we may estimate time delays. 
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